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Abstract 

OO ■ 

In this paper we study the asymptotic properties of the power variations of stochastic processes of 
the type X = Y + L, where L is an a-stable Levy process, and Y a perturbation which satisfies some 
I mild Lipschitz continuity assumptions. We establish local functional limit theorems for the power variation 

processes of X. In case X is a solution of a stochastic differential equation driven by L, these limit theorems 
provide estimators of the stability index a. They are applicable for instance to model fitting problems for 
paleo-climatic temperature time series taken from the Greenland ice core. 
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C3 



Stochastic differential equations are being used for quite a while as meso-scopic models for natural phenomena. 
In one of their simplest variants, they consist of deterministic differential equations perturbed by a noise term. 
The subclass in which the noise is Gaussian arises for instance from microscopic models described by coupled 
systems of deterministic differential equations on different time scales, in the limit of infinite fast scale, as the 
fluctuations of the slow scale component around its averaged version are considered. With a view in particular 
towards the mathematical interpretation of financial time series, the theory for stochastic differential equations 
the noise term of which is given by general (discontinuous and non-Gaussian) semimartingales has received 
considerable attention during the recent years. 

It is reasonable on a quite general level to model real data by stochastic differential equations. Usually 
. neither their deterministic nor their noise component can be deduced from first principles for instance from 
7-H ' microscopic models, but may be selected by statistical inference or model fit from the time series they are 
supposed to interpret. The central question of the model selection problem that motivated this paper asks for 
the best choice of the noise term. 

More formally, suppose we wish to model a real time series by the dynamics X — (X t )t>o of a real valued 
process of the type 



X t = x+ [ f(s,X s )ds + r) U t>0, (1.1) 
Jo 

where the process 77 is the noise component. Then the problem of a model fit consists in the choice of a drift 
term / and a noise term -q, so that the solution of (jl.l|) is in the best possible agreement with the data of the 
given time series. 

The example of (|1.1|) which inspired us most was investigated in the papers [21 0] , where P. Ditlevsen used 
(jl.ljl as a model fit for temperature data (yearly averages) obtained from the Greenland ice core describing 
aspects of the evolution of the Earth's climate during the last Ice Age, extending over about 100,000 years, 
in particular the catastrophic warmings and coolings in the Northern hemisphere, the so-called Dansgaard- 
Oeschger events [5] . The drift term was chosen as a gradient of a double well potential (climatic quasi-potential) , 
the local minima corresponding to the cold and warm meta-stable climate states. In order to find a good fit for 
the noise component, P. Ditlevsen performed a histogram analysis for the residuals of the ice core time series, 
the temperature increments measured between adjacent data points, i.e. years. He conjectured that the noise 
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may contain a strong a-stable component with a 1.75, and plotted an estimate for the drift term assuming 
the stationarity of the solution. 

In this paper we resume Ditlevsen's model selection problem for the fit of the noise component from the 
perspective of a new testing method to be developed. Following Ditlevsen, we work under the model assumption 
that the noise ij has an a-stable Levy component (symmetric or skewed). We search for a test statistics 
discriminating well between different a, and capable of testing for the right one. We shall show that this job is 
well done by the equidistant p-variation — in the sequel called p- variation — of the process X defined for all 
p > as 

[nt] 

V?[X) t := \ A i X \ P > 

2=1 

where A™X := X± — Xi-i for 1 < i < n, n > 1. We first observe that for those values of p relevant for our 
analysis the main contribution to the p-variation of X comes from the a-stable component of the noise. We next 
prove local limit theorems which hold under very mild assumptions on the drift term / and allow to determine 
the stability index a asymptotically. We finally use these limit theorems in Section [3] below to analyse the 
real data from the Greenland ice core with our methods, and come to an estimate a ~ 0.7, surprisingly quite 
different from the one obtained by Ditlevsen [3] 0] . 

The paper is organised as follows. In section [2] we set the stage for stating our main results, which are 
applied to the ice core data in section [31 The remaining sections are devoted to the proof of our main functional 
limit theorems, starting with the convergence of the finite dimensional laws in section H] continuing with the 
proof of the tightness of the laws in the Skorokhod topology in section [5] and ending with the robustness proof 
of the convergence with respect to adding terms of finite variation in section [5J 

In this paper, denotes convergence in the Skorokhod topology, '— »' denotes convergence of finite- 
dimensional (marginal) distributions, and ' stands for uniform convergence on compacts in probability. 
We denote the indicator function of a set A by Ia, and A denotes the complement of the set A. 



2 Object of study and main results 



Let (fi, T, {Tt)t>Qi P) be a filtered probability space. We assume that the filtration satisfies the usual hypotheses 
in the sense of [12], i.e. it is right continuous, and J-q contains all the P-null sets of T . 

For a E (0,2] let L = (L t ) t >o be an a-stable Levy process, i.e. a process with right continuous trajectories 
possessing left side limits (rcll) and stationary independent increments whose marginal laws satisfy 



lnEe iALf = 



-tC a \\\ a (l -i/3sgn(A) tan — 
-tC|A|(l-i/?-sgn(A)log|A|), 

7T 



a = 1, 



t > 0, 



(2.1) 



C > being the scale parameter and (3 E [—1,1] the skewness. We adopt the standard notation from [13] and 
write Li ~ S a (C,P,0). 

We also make use of the Levy-Khinchin formula for L which takes the following form in our case (see [5] 
Chapter XVII. 3] for details): 



InEe 



iALi 




1 - iAxI 
1 — iA sin x 



{|x|<l} 

1 



4a:<0} 



1 + /9 T 



Hx>0} 



dx 



\l+a ' 



•0T 



2 Hx<o} ■ 2 - 

where Cf denotes the scale parameter in Feller's notation and equals 

C Q (cos(^)r(a)) _1 , a^l, 
C-, a=l. 



{x>0} 



dx 



ll+a ' 



a^l, 
1, 



(2.2) 



a 



C F = 



(2.3) 



Recall that a totally asymmetric process L with /3 = l(/3=— 1) is called spectrally positive (negative). A 
spectrally positive a-stable process with a E (0, 1) has a.s. increasing sample paths and is called subordinator. 

The main results of this paper are presented in the following three theorems. The first theorem deals with 
the asymptotic behaviour of the p-variation for a stable Levy process itself. As we will see later, this behaviour 
does not change under perturbations by stochastic processes that satisfy some mild conditions. 
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Theorem 2.1 Let (L t ) t >o be an a-stable Levy process with L\ ~ S a (C, /3, 0). //p > a/2 i/ien 

(V(i) t -ntB n (a J p)) t > ^(Lj)t>o as n -> oo, (2.4) 
where L[ ~ S a i p (C' ', 1,0) wzi/i </ie sca/e 



C"= 4 



/COS ^ r 1 - ^) X p/a 

Vcos(^f)r(l-a)J ' (2.5) 
C, a = p. 

Tfte normalising sequence (B n (a 7 p)) n >i is deterministic and given by 

(n-^BlL^, pe(a/2,a), 
B n (a,p) = lEfan(n- 1 \L 1 \ a ), p = a, (2.6) 
[0, p > a. 

We remark that the skewness parameter (3 of L does not influence the convergence of V p n (L) t and does not 
appear in the limiting process since the p- variation depends only on the absolute values of the increments of L. 
Moreover, for p > a the limiting process L' is a subordinator. 

We next perturb L by some other process Y. We impose no restrictions on dependence properties of Y and 
L. The only conditions on Y concern the behaviour of its p- variation. We formulate two theorems, the first for 
p £ (a/2, 1) U (a, oo), and the second for p £ (1, a]. 

Theorem 2.2 Let (L t )t>o be an a-stable stochastic process, with L\ ~ S a {C, /3,0) and (Y t )t>o be another 
stochastic process that satisfies 

V; l (Y) u ^ p -0, n^oo, (2.7) 

for some p £ (a/2, 1) U (a, oo). Then 

(Vp(L + Y) t — nt B n (a,p)) t >o ^ {L' t )t>o as n — > oo, (2.8) 

with L' and (B n (a,p)) n >i defined in (|2.5[) and (|2.6D . 

The methods used to prove Theorem 12.21 do not work for p £ (l,a], and in this case we have to impose 
stronger conditions on the process Y. 

Theorem 2.3 Let (L t )t>o be an a-stable stochastic process with L\ ~ S a (C, /3, 0), a £ (1,2) and let (Y t )t>a be 
another stochastic process. Let p £ (l,a] and T > 0. If Y is such that for every 8 > there exists K{8) > 
that satisfies 

P(|y s (w) - Y t {uj)\ < K(6)\s-t\ for all s,t £ [0,T]) >l-5, (2.9) 
the process Y does not contribute to the limit of V™(L + Y), i.e. 

{Vf(L + Y) t -ntB n (a,p)) < t < T °{L' t ) < t < T , n - oo (2.10) 
wii/i L' and (_B„(a,p)) n >i defined in (|2.5p anc? (|2.6p . 



To be able to study models of the type p.l[) we formulate the following corollary of Theorems 12.31 and 
which takes into account that Lebesgue integral processes are absolutely continuous w.r.t. the time variable and 
thus qualify as small process perturbations in the sense of the Theorems. 

Corollary 2.1 Let (L t ) t >o be an a-stable stochastic process, with L\ ~ S a (C, /3,0), and let f : K + xI^M be 
a locally bounded function such that for some x £ R and T > the unique strong solution for 

X t =x+ f f(s,X s )ds + L t (2.11) 
Jo 

exists on the time interval [0, T]. Then for p > a/2 we have 

{V p n (X) t - ntB n (a,p)) <t<T ° (L' t ) < t < T (2.12) 
as n — > oo, with U and {B n (a,p)) n >\ defined in (|2.5[) and (|2 .6[) . 
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The functional convergence of power variations of symmetric stable Levy processes to stable processes was 
first studied by Greenwood in 6 , where more general non- equidistant power variations were considered, and 
in particular for p > a the convergence to subordinators was proved. Further, more general results on power 
variations of Levy processes are obtained by Greenwood and Fristedt in [7]. In [5] and [5], Jacod proves 
convergence results for p- variations of general Levy processes and semimartingales. In particular, several laws 
of large numbers and central limit theorems are established. Our results are different from Jacod's because we 
consider processes possessing no second moments so that only the generalised central limit theorem can apply. 
Moreover, we consider in addition convergence of perturbed processes. Corcuera, Nualart and Woerner in [T] 
consider p- variations of a (perturbed) integrated a-stable process of the type X = Y + f u s- dL s with some 
cadlag adapted process u. For u = 1, our setting results. The paper [Tj contains a law of large numbers for 
< p < a and a functional central limit theorem for < p < a/2. However, very restrictive conditions on 
possible perturbation processes Y are imposed, so that the results are not applicable to processes of the type 

(ED. 

3 Applications to real data 

In this section we illustrate our convergence results and show how they can be applied to the estimation of the 
stability index a. We emphasise that the conclusions obtained are somewhat heuristic. Additional work has to 
be done to provide more precise statistical properties of the p- variation processes as estimators for the stability 
index, and to describe the decision rule of the testing procedure along with its quality features. 

We first work with simulated data. Assume they are realizations of the SDE (|2. 1 1|) where L is a stable 
process with unknown stability index a. From the data set we extract m samples with n data points each by 
taking adjacent non-overlapping groups of n consecutive points. This way we get the samples (A"W)i<i< m with 
XW = (X± , , . . . , X^), 1 < i < m. Along each sample we calculate the p-variation 

n 

V^{X^) 1 =^2\A]X^\ P , l<i<m, 

3 = 1 

where the parameter p takes values in some appropriately chosen interval [p\ , p 2 ] ■ 

Due to Corollary the random variables F p "(XW)i converge to the stable random variable L[, possessing 
stability index a/ p. In order to estimate a, we compare the law of V^(X^)i with some known stable reference 
law. Thus for each p we calculate the empirical distribution function of (V//(X^)i)i<i< n given by 

1 m 

G p , n (x) ^-^If-oo,]^ 1 ")!), xel, (3.2) 
i=i 

and consider for convenience the reference ^-subordinator with scale parameter C > whose probability 
distribution function Fi/ 2 ,c can be explicitly calculated by 

[qF e -C'/2y 

w*) = v^/ — dy > (3 - 3) 

This will be a candidate for the limiting law of the random variable L[. The scale parameter C — C'(C) of L[ 
is connected with the scale parameter C of L by the relation (|2.5j) . To determine the unknown value of C we 
numerically minimise in p G [pi,P2], C 6 [Cij C2] the following distance of the Kolmogorov-Smirnov type: 

D n (C,p) = sup \G p , n {x) - F 1/2 ,c(C){x)\, 

x>0 

where < p\ < pi and < C\ < C'2 have to be chosen appropriately. Assume D n (C,p) attains its unique 
minimum at C = C* and p — p* . Since V™{X)i converges to a ^-subordinator if and only if p = 2a, 
we immediately obtain estimates for the scale and stability index of the driving process L, namely C* and 
a* =p*/2. 

To test this method, we simulate m — 200 samples of the data from equation (|2.11[) with f(-,x) — cosx, 
x G R, and L\ ~ 50.75(6.35,0,0), n = 200. We find that the Kolmogorov-Smirnov distance D n (C,p) attains a 
unique global minimum at C* « 6.35 and p* « 1.5 corresponding to the true values of a and C (see Fig. [1]). 

We next study the real ice-core data, analysed earlier by Ditlevsen in |3, 4j. The log-calcium signal covers 
the time period from approximately 90 150 to 10 150 years before present. We divide it into m = 282 samples, 
each containing n = 282 data points. Then the Kolmogorov-Smirnov distance is minimised numerically over p 
and C according to (|3.4j) . 

4 



(3.1) 



(3.4) 



0.0 0.5 1.0 1.5 

Figure 1: D n (C* ,p/2) for the simulated data, L is a 0.75-stable Levy process, n = m — 200 




Figure 2: D n (C* ,p/2) for the ice-core data set, n = m = 282. 




Figure 3: D n (C,p/2) for the ice-core data set, C = 3.28, n = m = 282. 
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It turns out that D n (C,p) for the real data also exhibits a unique global minimum in the (C, p)-domain, 
which yields the estimate a* k, 0.7 for C* « 7.2, D n (C*,p*) rj 0.1 (see Fig. [21). It is striking that our estimate 
differs from Ditlevsen's by a quantity very close to 1. This discrepancy can be explained as follows. It turns 
out that the function p t— » D n (C, p) has two local minima for some values of C different from the optimal value 
C*. For example, for C = 3.28 there are two local minima at a± « 1.02 and «2 ~ 1-76 with corresponding 
distances D n (C, 2a\) « 0.178 and D n (C, 2a-i) ~ 0.265. Unfortunately the paper [4] only contains the estimated 
value of the stability index a ~ 1.75 of the (symmetric) forcing L, and not its scale. It is possible that under 
some a priori assumptions on C, Ditlevsen's method provides a locally best fit which is not globally optimal 
(see Fig. 13]). 

4 Convergence of the finite dimensional laws of V p n (L) t 

To prove the convergence of the marginal distributions we use the following theorem which is a direct result 
of the well known generalised central limit theorem for i.i.d. random variables with infinite variance (e.g. see 
Theorem 3 in Feller Chapter XVII. 5]). 

Proposition 4.1 Let (r?i)i>i be a sequence of non- negative i.i.d. random variables with a regularly varying tail 
such that 

P(m > x) » C 2 —^-x- a l v asx-> +oo. (4.1) 
a 

for some a 6 (0, 2), p > a/2 and C > 0. Then for any t > we have 

(-) a J2 T )i-hn(a,p)±t p /°'Z, asn^w, (4.2) 



with 

(f \ p/a 
-) E? 7i> pe(a/2,a) 

4„Ksm(^!). P n. ( 4 -3) 



n 

0, p > a. 



where Z ~ S a / p (C, 1,0) with C as defined in (|2.5p . 

Let L be an a-stable Levy process as defined in (I2.1[) and let p > a/2. To study the finite dimensional 
distributions of V p n (L) t we note that due to the independence of increments of L it suffices to establish the 
convergence of marginal laws for a fixed t > 0. Further, the stationarity and independence of increments of L 
and the self-similarity property L t = t l t a L\ implies that 

[nt] [n*] lAlrip 



1 = 1 8=1 



with A}L = L% ~ S a (C,P,0) being i.i.d. random variables. It is easy to see that the random variables |Afi| p 
have a distribution function with a regularly varying tail, namely 

P(|AjL| p > x) =P(|ii| > x 1 ^) «C f — i-"/' asx^+oo, (4.5) 

a 

and thus we can apply Proposition 14. II to the sum (|4.4p . Taking into account that nt/[nt] -> 1 as n -> oo, we 
obtain convergence of the finite dimensional laws in Theorem 12.11 



5 Tightness of the laws 
5.1 Aldous' criterion 

To establish the tightness of the sequence (V™(L)) n >i is more complicated. Although the idea of the proof is 
based on an application of Aldous' criterion, the technical details depend strongly on the relationship between 
p and a. We first formulate a version of Aldous' criterion for tightness that is applicable in our case. 

On the probability space (fi, J 7 , P) we define nitrations Q n = (£/™)t>o generated by the power variation 
process Vp(L), i.e. 

Qt = °(V p n (L) t : s<t)= a{\L± | : < k < [nt]). (5.1) 



G 



Let T£j be the set of ^"-stopping times that are bounded by N > 0. Then the sequence V™(L) is tight if and 
only if the following two conditions hold (see p. 350 and p. 356 in [TO]): 

1. for all e > and N > there are no > 1 and K > 0, so that for all n > no we have 

P(supV p n (L) t >K)<e; (52) 

t<N y ' ' 

2. for all N > and e > we have 

limlimsup sup P(|F p "(L) T - V p n (L) s \ > e) = 0. 

010 n-»oo S\T£T™, (5.3) 

s<t<s+0 

Condition (|5.2p is obviously satisfied due to the monotonicity of the p- variation process and convergence of 
its marginal distributions to a stable law. 

To deal with condition (|5.3p we note that since (Vp(L) t )t>o has piecewise constant paths, the nitrations Q n 
increase only at discrete time instants t = k/n, k > 0, i.e. 

0« = 02 for *<t<*±± (5.4) 
n n 



Thus instead of Tj^ we can consider in (|5.3p a family of ^"-stopping times with finitely many values 

5ft :={SeK:S takes values in {0, -, . . . , N}}. (5.5) 

n 

5.2 Tightness for p > a 

The case p > a is simple because the compensating sequence B n j(a,p) vanishes. Thus we can use the mono- 
tonicity of the process V™(L) t . 

Let N > 1 and e > be fixed. We use the stationarity and independence of increments of L and convergence 
the marginal distributions of V™(L) to L' to obtain the following limit: 

limlimsup sup pf \V"(L) T - V?{L)s\ > e] 

010 n->oo S,TeS%, ^ ' ' 

S<T<S+0 

nT 

= limlimsup sup p( V |A"L| P > e) 

810 n ^oo S ,T65S. i=t<?+l 

s<T<s+e 

nS+[n9] 

< limlimsup sup pf V \A?L\ P > e) 

N 2—nS+l 
nN 3+[n6] 

= limlimsup sup Vp( V I A"L\ P > e, S = - 
= limlimsupPI V |A"L| P > e) 

2 — 1 

= limlimsupP(^"(L) e > e) 

010 n^oo 

= limP(L' e > e) = 0. 

010 

5.3 Tightness for ~ < p < a 

The case p < a has to be treated differently since the process 

X" = V p n (L) t - [niln-v'^L^ (5.7) 

need not be monotone. 

Let 6 > 0. As in the preceding section we use the independence and stationarity of increments of L to obtain 



(5.6) 
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the estimate 



nS+k 



sup P(|X£ > e) < sup Pf max I V (|A"L| P 
<;» Se5™, ^ l<fc<[nfll I ^— ' V 



S<T<S 



-n- p / a E|Li| p 



•i_-nS'-t-l 



> e 



<P( max |Vf|A"L| p 



,-pA 



E|Li| p 



< max 3P( \y (|A"L| 



E|L!| P 



> e 



> 



(5.8) 



37' 



where the last inequality follows from Lemma 20.2 in [13] . Now we have to show that the probabilities in the 
latter inequality converge to zero uniformly in k as n — > oo and 010. 
For any 8 > we can choose 9q > and n > 1 such that 



[nfl 



'(I X (i A r^r -rr p/a E|Lii p 



i=l 



> - < - 

" 6/ - 6 



(5.9) 



for all n > ^f. Without loss of generality we can assume that 6*n < 1- We will show now that 



•M ^ MAf£| p - n~ p ^ a E\Li\ p 



> - < - 



(5.10) 



holds for all 9 < 9o and n > uq. Indeed let us fix 9 < 6q and assume first that 9 < 9q/2. With the help of the 
triangle inequality and the estimate (|5.9p we can conclude that for n > uq 



[nfl] 



(|A£\L| p -n- p / Q! E|£il p 



> 



[nfl Q ] 



< P 



4 + * 



(|£ (ML\*-n-*' a E\L l \* 

[n6 ]-[n6] 

J2 {\KL\ P -n- p/a nLi\ P )\>l) 



9 ] 

> 

i— l+[n- 



(| Y, (|A?£| p -rr p/a E|£ 1 l p 



£ 

> - 

~ 6 



(5.11) 



Recalling that 9o — 9 > 6q/2 and 9a < 1 we find n > n' > % such that [n'flo] = [ n 9o] — [n9]. Thus employing 
the self-similarity of stable processes and (I5.9[) we have 



[nfl ]-[nfl] 



e 

> - 
- 6 



'(J ]T (|A?£|*-n-»/°'E|L 1 p 

i=l 
[n'flo] 

y (|A"L| P -n- p / Q E|Li| p 

/ / [n'flo] 

' P/ "| £ (|A?'L|" - (nT^ElLx 



e 

> - 
~ 6 



(5.12) 



> 



[n'flo] 



<p(| (\ A i' L \ P - (n')" p/a E|£i| p 



> - < -. 

" 6/ - 6 



So the inequality (|5.10jl is proved for fl < flo/2 and n big enough. Further, if flo/2 < 9 < 9q the proof is even 
easier because for n > no we can directly find n' with the same properties as above such that 



'(|^(JA?L| P -n-^ElL-V 



> 



Y (|A"L| P -n^ p / Q E|L!| p 



e 

> - 
~ 3 



[n'flo] 



[n'flo] 



(5.13) 



<P(| Y (M' L\* ~ (nT p/a nL 1 \*)\ >|) <P(| (|A?'L|* - (nT^E^ 



> T < 



8 



Thus we have demonstrated that for any S > and e > there exist # > and no > 1 such that 
3P(| XI (|A"L|f - n^ElLil^) |>|)<*, « > «o, B < B , 



(5.14) 



which together with Ij5.8|> implies the second condition of Aldous' criterion, namely 



limlimsup sup P(|X£ - X%\ > e) < (5. 

s<T<s+e 



(5.15) 



But <5 was arbitrary. 



5.4 Tightness for p = a 

For this part of the proof we need two lemmas that will enable us to control moments of stable processes. 

Lemma 5.1 Let L be an ot-stable Levy process and let k > 1. Then there are positive numbers C a ^, K a ^ and 
N a such that for y > N a we have 



K a . k y^ a < EdLxl^I^^}) < C a , k y^ a 

if k > 1 and 

K a ,i lny< Ediil^ii,!^}) < C a<1 Iny. 
In particular there is C a > such that P(|Li| > y) < C a y~ a for y > N a . 



(5.16) 
(5.17) 



Proof: We first prove the estimates from above. Let F denote the probability distribution function of L\. 
Then there are positive constants c a and n a such that for x > n a 



1 - F(x) < and F(-x) < — 



(5.18) 



Integration by parts yields for y > n a 



E 



(\Li\ ka h\ Ll \<y}) = f \x\ ka dF{x) 



-x ka F(-x) ' + I kax ka - l F(-x)dx 
o ' 



■x ka (l-F(x)) + / kax ka - L (l- F(x))dx 
o 



< 2 / kax^- 1 dx + 2 I kax Ka ~ L ^ dx 

JO Jn a xDt 

= 2n k a a + 2c a f V kax^- 1 ^- 1 dx 



ka-1 L a 



(5.19) 



2n k a a + 2c a 



< 2n k a a + 2c a 



fc(lny — lnn a ), k = 1, 

'_A_ y (fc-l)a ; fc>1) 

klny, k = 1. 



Now we can choose constants C Qj fc and N a such that the estimates from above are satisfied. The estimates from 
below follow analogously from the inequalities 



k k 
1 - F(x) > — and F(-x) > — 



(5.20) 



for some k a > and x big enough. 



Lemma 5.2 Let L be an a-stable Levy process. For any a,b > there exist C a ^ n (a,b) and N a (a,b) such that 
for A 6 [0,2/a] and n > N a (a,b) we have 



EexpMA^ 



exp (iAEsin (^-)) <C a . n (a,b) 



with C a . n (a,b) satisfying 



lim lim C a ^ n (a, b)nb = and lim lim C a n (a, b)n = 0. 

b— >0 n—>oc 



a — >oo n — >oo 



(5.21) 



(5.22) 
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Proof: We split the left-hand side of (|5.21[) into the real and the imaginary part to obtain the simple estimate 



Eexp (iA^) - exp (iAEsin (^)) 



< 



Ecos 



( MLi\ a 
\ n 



— cos ( AE sin 



\Li 



E sin -sin (\Vsm l -±- 

\ n ) V n 



\Li 



(5.23) 



(5.24) 



Throughout the proof we shall use the following elementary inequalities: 

x 3 . 

x — — < Sim < x, x > 0, 
x 2 

cos x > 1 — — , x € M. 
Let a, b > 0, A S [0, 2/a] and denote m n = m n {a,b) := na 1 / 4 ?) 1 / 2 . 

We first estimate the real part. To this end, we apply Lemma 15.11 with n big enough so that ml/ a > AT a or 
n > A^a" 1 / 4 ^ 1 / 2 , to obtain 

l>Ecos(^L!) 



> Ecos 



>E 1 



A 2 |Li 



= 1 -E 



2n 2 
A 2 |Li| 2 " 



L {\L 1 \< m L n ' a } 



P(|ii|>my«) 



^Vi<^a)- 2P(|Ll| ^ m » /a) 



> 1 



> 1 



C Qj 2A 2 m rl 2C Q 
2n 2 m n 
2C Q , 2 o 1 / 2 



2C a 



na 7 / 4 na 1 /^ 1 / 2 
Analogously we get an estimate for the second summand: 

'|£i| a \\ . , A 2 



l-CW(o,6). 



l>cos(AEsin(^)) >l-^-( E sin( 



> 1 - 



A 2 
2 

A 2 



Esin 



ll {|ii|<mJ > /a! } 



>1- T (E 



Summarising we have 



> 1 
= 1 

Ecos 



n {|iil<< /Q } 
2 /C a iln 



\Lx\ a 
n 

P(|Li|>my°i 
P(|ii|>my«) 



4(- 

^(C^lnina 1 '^ 1 ' 2 ) , C Q 
a 2 V 



AEf 



na 1 /46i/2 



l-C«(o,6). 



<CW(a,6) + cW(a,6). 



We now estimate the imaginary part. With analogous arguments we obtain 



"(^)--(»-(*f)) 



l^il"' 



<»-.(^)-«- »( M f)^( E -( M f)) 



l^i| Q \ , A 3 



< E 



(( sin i^T-) - Asin (^K^i^i) + C 1 + A ) p (" Ll " * m ^ 

4 /Ca^ln^a 1 / 4 ^ 1 / 2 ) 
V na 



n 

a 



3a 3 



6 1 / 2 ) 



-A|Li| Q A|Li| Q , A|ii 



(( n n ()// :! 



3t> 



4 / Caaln^a 1 / 4 ^ 1 / 2 ) 
3a 3 V na na 1 / 4 ^ 1 / 2 / 



{|il|<m„ /Q } 
3 



(1 + 2/a)C a 



in,, 



V 



< 



C a , 3 6 (l + 2/a)C a 4 / Caqln^a 1 /^ 1 / 2 ) C a \3 (3) 

3noVa + naW 3a 3 V na + naW/aJ ' ° Q '" la ' j ' 



(5.25) 



(5.26) 



(5.27) 



(5.28) 
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The lower bound is obtained similarly, and reads 



Esin 



\\Li\° 



> Esin (^)-AE sin (£f) 



sin ( AE sin 



\Li\ a 



> E sin 



n 

MLi\ a 



A sin 



\Li\ 



> E 



(( 



n n 

la \3|r 13a 



{\Li\<mi /a } 



X\Lt\ a A 3 |L!| 3Q A|Li|° 



> 



n 6n 3 n 

4C a , 3 b (1 + 2/a)C a 
'3a 5 / 2 n na 1 /^ 1 / 2 



{\L 1 \<rn 1 J c "} 



-(l + X)P(\L 1 \>m 1 J a ) 
(1 + 1)^ 



(5.29) 



Combining these estimates with (|5.27l) and denoting C a . n (a, b) := Ca} n (a, b) + Ca}i(a, b) + Ca}n(a, b) + C^n(a, b) 
we obtain inequality (|5.2ip . Property (|5.22[) of the limits is straightforward. ■ 
Now we can show tightness for p = a. We use the inequality analogous to (|5.8[) which states that for any 
e > 

\Lj\ 
n 



i=nS+l 



»P, P (| E (|A„L 4 r-Esin(^))|>e) 



s,teS' 

S<T<S+ 



< max 3P 

Kk<\n6\ 



*Li \ a - Esin 



e 

> - 
~ 3 



Therefore we have to show that 
with 

X 



limlimsup max PflX^I > e) = 

810 n^oo l<k<[n$] 



1=1 



A„L,;| Q -Esir 



liil 



, fc, n > 1. 



The truncation inequality (see e.g. Theorem 5.1 in |llj ) provides an upper bound for the tail of X^: 
P(\X*\ > e) < | / ' ( 1 - Ee iAX ") dX < 2 sup f 1 - sREe iAX ") , 



-2/e 



0<\<2/s 



(5.30) 

(5.31) 
(5.32) 

(5.33) 



where Sftc denotes the real part of a complex number c. Since X% is a sum of i.i.d. random variables, its 
characteristic function can be factorised and we get for A > 



Ee 



Eexp ^iA^ 
exp 



\Li\ a . /|ii| a ' 
■ — ■ E sin f 



fclnEexp(iA(^^ 



Esin 



Denote 



exp 



exp 



fcln 



Eexp(iA^) 



(5.34) 



fcln 1 



exp(iAEsin(^)) 

Eexp(iA^) - exp(iAEsin(^)) 



exp(iAEsin(^)) 
Eexp(iA^) - exp(iAEsin(^)) 



exp(iAEsin(^)) 
and note that due to Lemma 15.21 for any e, 9 > we can estimate 



|V„,a| = |Eexp (iA^^) - exp (iAEsin (j^-)) \ < C a , n (e,l 



(5.35) 



(5.36) 



with lim e ^ hm„_ ) . 00 C a , n (e, 9)9n = 0. 

Let 9 > be fixed and let no(a,e, 8) be such that \ipn,\\ < 1 for ft > no(c\!,£, 0). Recalling the elementary 
approximations ln(l + z) = z + p(z) and e z = 1 + z + q(z) with [p(z)|, \q(z)\ < \z\ 2 , \z\ < 1, and the estimate 
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(I5.33P we get 

P(|X*|> £ )<2 sup (l-3?(exp[fcln(l+V„,A)])) 

Ae[0,2/e] 



= 2 sup ( - M fe^n.A + kp(lp n ,\) + q{klp n .\ + kp(lpn,\))) 
AG[0 : 2/e] V V 7 

< 2 sup ( k\t(>n,x\ + k\p(ip n ,\)\ + \q(kif} n ,\ + kp(tp n ^)) 

lcln iw \ 



(5.37) 



Ae[0,2/e] 

For n > max(no, N a (e, 6)) with N a (e, 9) as defined in Lemma IS~2l we already know that 

£#„, A | + fe|p(V>„,A)l < + l)l^n,A| + *(« + l)hM a < 2^(n + l)C a , n (e, 0). (5.38) 
Applying Lemma T5. 21 we conclude that for < 9o, n > max(no, N a (e, Oq)) and any e > 

lim limsup max P(|X*| > e) = 0, (5.39) 

0— *0 ri—>oo Kk< \n6] V ' 



so that condition (|5.3p of Aldous' criterion holds. 

To prove condition (|5.2p we proceed in exactly the same way as for (|5.3p , taking advantage of the fact that 
for fixed N > 1 the equation 

lim lira C a n (-^> N)Nn = f5 40) 

holds due to Lemma [57^1 

6 Generalisation to sums of processes 

We finally discuss the situation of Theorem 12.31 where besides the Levy process L another process Y is given. 
To see that V p {L) and V™(L + Y) have equivalent asymptotic behaviour we apply Lemma VI. 3. 31 in [TU] . 

Under the conditions of Theorems and it is enough to show that V^(L +Y) - V p {L) u ^> p ' as n — > oo. 
6.1 Equivalence for p < 1 

Let L be an a-stable Levy process, p G (a/2, 1] and let Y be such that for V p (Y) U ^ P ' 0, as n — > oo. Note that 

due to the monotonicity properties of V^"(Y), the latter convergence condition is equivalent to V^{Y)n — > 
for any iV > 1. Then a simple application of the triangle inequality yields the proof. In fact, for any TV > 1 we 
have 

nN 

sup \V p n (L + Y) t -V p n (L) t \ <]T||Ar(£ + y)p , -|ArL| p | 

0<t<N , 

1=1 (6.1) 

nN y ' 

< Y, |A?(L + y) - A"i| p = V p n (Y) N £ 0, n - oo. 
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6.2 Equivalence for p > a 

p 

Assume again that V p n (Y) t — > 0, n — > cx), for any t > 0. Denote m := [p], so that p = m + q with q G [0, 1). 
Then for any iV > 1 we have 

[nt] 



sup |y p "(L + Y) t - V p (L) t \ = sup | + ^)| m+9 - |A n L,p+« 

[nt] 



t<N 



m — k 



= su P |£ |a"(l +ni 9 EJ l A «^l"|A?y 

^E UjEH^^ + ^l 9 !^!"!^ 

fe=0 ^ ' i=l 
nJV 

+ 53 ||A?(L + y)|9|A?L| m - |A"L| m+ « 



m — k 



\A?L\ m+q 



i=l 
m — 1 



k=0 
nN 



nN 



m—1 / x nJV 

m 



< E ( r E |Ari| fe+9 |A:[ l rr i - fc + E U E |A?L| fc |A?yr- fe+ <? 



i—k 



+ J2 \ A ?L\ m |A?(L + Y)| 9 - |a?z| 9 

t=l 

m — 1 / \ TiiV 

^E (T Ei A ^i fc+9 i^ 

fc=0 ^ ' i=l 
nN 

+ I {q>0} Y / \^LH^Y\ q . 



k=0 



m — 1 



fe=0 



nN 



E : El^l fe |A«r 



rn— 



(6.2) 



The right-hand side of the latter inequality is essentially a sum of 2m+l terms of the type Y^i=i I A™£| fe | A™y \ p ~ k , 
where /c G [0,p] nN. 

Applying Holder's inequality we get 



it] 



it] 



53 iA„L 4 | fc |A„y 4 r fe < (53 |A n x i r) fe/p (53 



["*] 



i=l 



i=l 



(p-k)/p 



(6.3) 



The first factor in the latter formula converges in probability to a finite limit, since p > a. The second factor 
converges to in probability due to the assumption on Y, and Theorem 12.21 is proven. 



6.3 Equivalence for a G (1,2), p G (I, a] 

The main technical difficulty of this case arises from the fact that the p- variation of L for p < a does not exist. 
In particular, the events when increments of the stable process L become very large have to be considered 
carefully. For T > and some c > which will be specified later define the following sets: 

J?H := {i G [0, [nT]) : |A?L( W )| > c}, 

A n c (j) := G n : | J™ M| = j}, j = 0, . .., [nT]. 1 " j 

The set J™ contains the time instants (in the scale i), where the increments of the process L are 'large', i.e. 
exceed c. The set -A™ (j) describes the event that the number of large increments equals j. 

Let 5 > 0, e > 0. According to the conditions of Theorem (|2.3p . let A"' := AT(|) > 0, and set 

B = {uj: \Y B (u)-Y t (u))\ <K'\s-t\ for all M e [0,T]}, (6.5) 

so that P(B) > 1 - §. 
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We estimate 



( sup \V p n (L + Y) t -V p n (L) t \>e) 



iT] 



(6.6) 



<p({^||A?(L + y)|P-|A?L|p| >e}ns) + P(B) 

i=l 

[nTI [ ™ T] <5 

^ E p ({ E | i A "( L + y )i p - i A " L H > e } n A c (•?') n B ) + 2 

j=0 i=l 
[nT] 

^ E p ({ E |l A "( x + - l A " L l P | > | } n A " (i) n b) 

j=0 igj? 
[nT] 

+ E p ({ E \\ A ^ L + Y )\ P - \^ L \ P \ > \) n A c 0") n b) + g 

= :Z?( 1 )(n,c)+ J D( 2 )(n, C ) + ^. 

In the following two steps we show that for appropriately chosen c = c(e) > and n big enough, D^\n, c) = 
and D^ 2 \n, c) < 5/2. This will finish the proof. 

Step 1. To estimate D^\n, c), let lu e B. Using the elementary inequality \a p — b p \ < max{a, b} p ~ 1 p\a — b\ 
which holds for a, b > and p > 1 we estimate 



£ ||A?(L(o;)+y(a;))r-|A?L(a;)r|< £ p( c +^) P ~' 

< E p(c + ^)'-|A?yi 

< E K c+ -) 

i0J™(ai) 



A?(LM+y( W ))|-|A?L(a;)| 



(6.7) 



n / n 



where we have used that the 'small' increments of L are bounded by c. So we have 

[nT] 



for all n>ni with 



D«(n,c) = E P ({| E l A " X + A " y l P - l A " L l P | > |} n n S ) 

s£ p (M c + t)'' ^>|}n^(i)ns) 

s ;; Kc+ ^ r vr>i).o 



(6.8) 



ni := 



2^ 



and c = c(e) := if — ^— — 
v 7 A\2pK'T 



Step 2. To estimate _D^ 2 '(n,c), let n > ^- so that for u> G £? and i G </™(w) we have 



if' 

< — < 1. 

nc 



(6.9) 



(6.10) 



A?L(w) 

By means of the elementary inequality (1 + \x\) p — 1 < 3|x| which holds for 1 < p < 2 and \x\ < 1 this implies 
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that 



£ \M(L(u)+Y( U ))P-ML(u)P 

i£J?(ci>) 



< ^ |A?L( W )|f3 



A?i(«) 



3if 



(6.11) 



This in turn immediately yields 

DW(n,c) < n { £ M L \ P ^ > f }) * E P (^(i) H (J > ^-}). (6.12) 



i=i 



Since all A"L, z = 1, . . . ,[nT], are i.i.d., and only j of them exceed the threshold c, we can estimate the 
probability for this event precisely. Indeed, denoting 



Pn := P(|A?L| > c) = P(|Lx| > m 1/Q ), 
we continue the estimates in (|6.12|) to get 

IT) 

i 1 j 



(6.13) 



£ (2) (™, c) < £ ( [ri J ] )p({{|A?L| > c,i = 1, . . . , j} n {|A?L| < c,i = j + 1, . . . >T]}} n|j{|A^> ^ 7 
jP({|A?i| > c,t = 1, . . . ,i} n {|A"L < c,i = j + 1, . . . , [nT]} n {|A?L|f > ^}) 



[nT] 

= E 

[nT] 

^E 

3=1 



nT] 

i 



nT] 
j 



With help of the inequalities (|5.18|) . we obtain the estimates 



Pn < — 



and 



p(|A?£|> > ^ = Pf|ix| > f^V'VM < 



-a/p 



(6.14) 
(6.15) 
(6.16) 



holding for some c > 0, for all n > n% > — and 1 < j < [nT]. Combining (|6.14| and (|6. 16|1 . denoting the 
constant pre-factor by C, and recalling that ~ < 2, we obtain for n>rii that 

< ^7FE ( [r ; T] )i 1+Q/ M(i -aO^"' < ^rWE ( [n 7 ] )iV„(i- P n)^-'. (6.17) 

The sum in the previous formula represents the third moment of a binomial distribution, and thus can be 
calculated explicitly. By means of the asymptotic inequality (|5.18|) we get 



D {2) (n,c) < C 



(nT - 2)(nT - l)nTp% + 3(nT - l)nTp 2 n + nTp n 
p n n 1+a /P 



^ ^((nT Pn ) 2 + 3nTp n + l). 



(6.18) 



Now choose n > n$ > 712 big enough to ensure that this expression is smaller than |. This completes the proof 
of Theorem P 
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